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A posteriori error estimates for mixed finite element Galerkin 
approximations of second order linear hyperbolic equations 
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Abstract 

In this article, a posteriori error analysis for mixed finite element Galerkin approxi¬ 
mations of second order linear hyperbolic equations is discussed. Based on mixed elliptic 
reconstructions and an integration tool, which is a variation of Baker’s technique intro¬ 
duced earlier by G. Baker (SIAM J. Numer. Anal., 13 (1976), 564-576) in the context 
of a priori estimates for a second order wave equation, a posteriori error estimates of 
the displacement in L°°(L^)-norm for the semidiscrete scheme are derived. Finally, a 
first order implicit-in-time discrete scheme is analyzed and a posteriori error estimators 
are established. 


Key words, second order linear wave equation, mixed finite element methods, mixed 
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1 Introduction 

In this paper, we discuss a posteriori error estimates for mixed finite element Galerkin 
approximations to the following class of second order linear hyperbolic problems: 

(1.1) utt-V ■ (AVu) = f in O X (0,r], 

(1.2) = 0 u\t=o = uo and ut\t=o = ui- 

Here, Q C SA is a bounded polygonal domain with boundary dtl, 0 < T < oo, ut = ^ 
and A(x) = (aij(x))i<jj <2 is a symmetric and uniformly positive definite matrix. All 
the coefficients a^-’s are smooth functions of x with uniformly bounded derivatives in H. 
Moreover, the initial functions uq = uo{x), ui = ui{x) and the forcing function / = f{x,t) 
are assumed to be smooth functions in their respective domains. 

In recent years, there has been a growing demand for designing reliable and efficient 
space-time algorithms for the numerical computation of time dependent partial differential 
equations. Most of these algorithms are based on a posteriori error estimators, which provide 
appropriate tools for adaptive mesh refinements. For elliptic boundary value problems, a 
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posteriori error estimates are well developed (see, [31 [32]). Adaptivity with a posteriori 
error control for parabolic problems has also been an active research area for the last two 
decades (cf. [HI |M1 E511331 EIIHIIS] and references, therein). For the time discretization, 
some results on a posteriori error estimations for abstract first order evolution problems are 
available in the literature (cf. [H EH EEl ESI [5U] 1. 

In the context of second order wave equations, only few results are available on a pos¬ 
teriori error analysis, see, PI H [1311131171131]. Further, it is observed that the design and 
implementation of adaptive algorithms for these equations based on rigorous a posteriori 
error estimators are less complete compared to elliptic and parabolic equations. Based on a 
space-time finite element discretization with basis functions being continuous in space and 
discontinuous in time, a priori and a posteriori error estimates for second order linear wave 
equations are proved in pi]- Asymptotically exact a posteriori estimates for the standard 
finite element method are proposed and analyzed in mm by solving a set of local elliptic 
problems. The recent results in [T] EQ] cover only first order time discrete schemes. In [7], 
the second order wave equation is written as a first order system and a first order implicit 
backward Euler scheme in time is used with continuous piecewise affine finite elements in 
space. Further, rigorous a posteriori bounds have been established using energy arguments 
and adaptive algorithms based on the a posteriori bounds are discussed. In pO], based on 
Baker’s technique a posteriori bounds are derived for the semidiscrete in L°°(L^)-norm and 
for first order implicit-in-time fully discrete schemes in i°° {L‘^)-noTm. The fully discrete 
analysis relies crucially on a novel time reconstruction satisfying a local vanishing-moment 
property, and on a space reconstruction technique used earlier in p8] for parabolic problems. 
In [13], an adaptive algorithm in space and time which is based on Galerkin space-time dis¬ 
cretizations leading to Newmark scheme is analyzed. Further, goal oriented a posteriori 
error estimates are derived and some numerical results are provided to demonstrate the 
efficiency of error estimators. In [3T], the author has studied an anisotropic a posteriori 
error estimate for a finite element discretization of a two dimensional wave equation. The 
estimate is derived in the T^(0, T, 77^(n))-norm and it turns out to be sharp on anisotropic 
meshes. 

For higher order time reconstruction for abstract second order evolution equations, one 
may refer to the recent papers [23] [22]. In p3], an adaptive time stepping Galerkin method 
is analyzed for second order evolution problems. Based on the energy approach and the 
duality argument, optimal order a posteriori error estimates and a posteriori nodal super¬ 
convergence results have been derived. An adaptive time stepping strategy is discussed 
and some numerical experiments are conducted to assess the effectiveness of the proposed 
scheme. In a recent work [22], second order explicit and implicit two-step time discretiza¬ 
tion schemes such as leap-frog and cosine methods are discussed and a posteriori estimates 
using a novel time reconstruction are derived. Further, some numerical experiments are 
conducted to confirm their theoretical findings. 

For space-time adaptivity, the finite element discretization depends on the space-time 
variational formulation and its error indicators include both space and time errors. Recently, 
attempts have been made to exploit elliptic reconstruction to prove optimal a posteriori 
error estimates in finite element methods for parabolic problems |28j . In fact, the role 
of the elliptic reconstruction operator in a posteriori estimates is quite similar to the 
role played by elliptic projection introduced earlier by Wheeler [M] for recovering optimal 
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a priori error estimates of finite element Galerkin approximations to parabolic problems. 
This analysis is, further, developed for completely discrete scheme based on backward Euler 
method [2B], for maximum norm estimates m and for discontinuous Galerkin methods for 
parabolic problems [21]. In recent works [29] and m, the analysis is further extended to 
mixed FE Galerkin methods applied to parabolic problems. 

In this article, an a posteriori analysis is discussed for mixed hnite element Galerkin 
approximations of a class of second order linear hyperbolic problems. One notable ad¬ 
vantage of mixed finite element scheme is that it offers a simultaneous approximations of 
displacements and stresses, resulting in better convergences rates for the stress variable. 
This property is important in applications such as in the modeling boundary controllability 
of the wave equation. In the first part of this article, a semidiscrete scheme is derived us¬ 
ing mixed finite element method in spatial direction, while keeping time variable constant. 
Based on mixed elliptic reconstructions presented in [29], which depend explicitly on resid¬ 
uals and a time integration tool, a variant of Baker’s technique, a posteriori error estimates 
in L°°(L^)-norm are derived for the displacement u. For the time discretization, the time 
discrete scheme with the time-reconstruction proposed in [20] is applied and then, using 
summation tool, a posteriori error estimators in t'°°(L^)-norms are developed. Compared 
to m, our analysis is not only for mixed finite element method, but also it differs from 
the analysis of [20] in the sense that a time integration tool is used for deriving a 

posteriori estimators, as against the time testing procedure of Baker [6] used in m- 

The outline of this article is as follows. In Section [21 we introduce both weak primal and 
mixed formulations for the hyperbolic problem (ll.ip - (jl.2p and establish their equivalence. 
Section [3| deals with mixed elliptic reconstruction techniques proposed in [29] and a poste¬ 
riori estimates for the semidiscrete problem for both displacement u and its stress cr in 
L°°(L^)-norms are derived. Based on a first order backward differencing implicit method, 
a completely discrete scheme is proposed and related a posteriori error estimators are es¬ 
tablished in Sectional Finally, results are summarized in Section [5] with a brief outline on 
future work. 

2 On primal and mixed formulations 

We use the usual notations for the L^,Hq and spaces and their norms and semi¬ 
norms. Let be the dual space of Hq and let (•, •) be a duality paring between and 
Hq. Since we shall be dealing with time-space domain, we further introduce for a Banach 
space X with norm || • ||x, the space LP(0, T ; X) denoted by ]X{X), for 1 < p < oo with norm 
II ■ \\lp{X)- Moreover, denote by H'^{X) the space of vector valued functions cj) '■ (0, T) —^ X 

such that € L‘^{X) for j = 0,1,... , m with the standard norm || • \\h'^{X)- 

For the weak primal formulation, define a bilinear form for w, z € Hq 

a{w,z) := {AXw,Xz). 

Given / G uq G Hq and ui G L^, the weak formulation of (ll.ip - (ll.2p is to seek 

a function u : (0, T] —>■ Hq with m(0) = uq and ut{0) = ui such that 

(2.1) {utt,w)a{u,w) = {f,w) 'iw€HQ. 
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Note that for / E L^(L^), uq E Hq and ui E L^, there exists a unique weak solution u of 
(HHi-dril) satisfying u E L‘^{HQ),ut E L^(L^) and uu E Lp‘{H ^). Moreover, the equation 
m is satisfied for almost all t E (0,T]. For a proof, refer to Evans (|IS], PP- 399-408). 
For mixed formulation, let 

H(diu, n) = {(t)e {L^{VL)f : V • (/) E L‘^{n)} 

be a Hilbert space equipped with norm ||0||H(dii;,o) = (ll^lP + l|V ' 0lP)^- 
Now, introduce 


(2.2) <T = —AVu, 
and a = A~^. Then, the equation (II.ip is rewritten as 

(2.3) a<T-|-Vu = 0, utt + '^-cr = f. 

Set W = LP‘{VL) and V = H((izu,n). For given / E L^(bF), E VF, a weak mixed 

formulation for ()l.ip - ()1.2p is to find {u,cr) : {0,T] ^ W xV with u(0) = uq and ut(0) = ug 
such that (u, <t) satisfies 


and 

u,uteL\W), uueL\H-^) 

and a E T^(V) 

(2.4) 

{acr, v) — (ti, V • v) = 0 

VvE V, 

(2.5) 

{utt,w) + (V ■ a,w) = {f,w) 

Vw€W. 


Below, we discuss the equivalence of weak primal and weak mixed formulations. 

Theorem 2.1. The pair {u, a) E L‘^{W) x L^(V) with uu E and (uq, ui) E IF x VF 

is a solution of the mixed formulation (1231)-dZIl) if and only if u is a solution of the weak 
formulation dm) and a = —AS/u with u E L'^{Hq) and uq E Hq. 

Proof. Let (tt,cr) E L‘^{W) x T^(V) with uu E be a solution of ()2.4p - (|2.5p and let 

(p E T>(f4). Choose v = Curl cp '■= {—d(p/dx 2 ,d(p/dxi) in (|2.4I) . Since it is divergence free, 
we obtain 

(acr, Curl </>) = 0 a.e. t E (0, T). 

Using distributional derivative, it follows that 

(Curl (acr), cp) = 0\/(p € ^(ff), a.e. t E (0, T), 


and hence. 


Curl (A ^cr) = 0 in VfQ), a.e. tE(0,r). 
Now, a use of Helmholtz decomposition yields for some ip E Hq 

A~^a = Vip a.e. t E (0,T). 
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Apply this in (12.4p to arrive at 


(2.6) (V ■ V, V’ + rt) = 0 Vv G V, a.e. t € (0, T), 

which shows u = — G Hq. Hence for a.e. t G (0,T), 

A~^a = —Vu 


and 

<T = —AVu a.e. t G (0,T). 

On substitution in (|2.5I) yields (|l.ll) and hence, it satisfies (|2.ip for a.e. t G (0,T). Since, 
u,ut G LP‘{W), u G r] and (|2.6I) holds for t = 0. Thus, uq = —V'(O) G Hq. 

For the converse, let u be a weak solution of o satisfying dlT]). Now, set cr = 
—AVu G Then, the equation (II. ip becomes 

utt + V-cr = f in V'{Q.), a.e. tG(0,T). 

Since, / G L^(L^) and V • cr G therefore, + V • cr G and (|2.5I) is satisfied 

for a.e. f G (0,T). Note that acr = —A“^cr = —Vu. Multiply by v G V, and integrate over 
Q, to show that (|2.4p is satisfied. This concludes the rest of the proof. □ 

Since the weak primal formulation (j2.1h is well-posed, by the Theorem l2.1l on equivalence, 
the weak mixed formulation is well-posed. As a byproduct, since / G L?‘{LP‘) and V • cr G 
therefore, uu G 

Given A uniformly positive definite, there exist two positive constants oq and oi such 
that 

ao||<^|| < ll'^IU-i < oi||<^ll) where ||cr||^_i := (acr,cr). 

3 A posteriori error estimates for the semidiscrete scheme 

This section focuses on a mixed finite element method for the hyperbolic problem (II.ip - 
(HH) and a posteriori error estimates are derived for the semidiscrete mixed Galerkin 
approximation to (ll.ip - (11.2h . 

For the semidiscrete mixed formulation corresponding to (I2.4I) - ()2.5I1 . let Th = {AT} be a 
shape-regular partition of the domain O into triangles of diameter hx = diam(iF). To each 
triangulation Th, we now associate a positive piecewise constant function h{x) defined on 
O by /i|/c = hx VAT G Th- Let F/j denote the set of all internal edges E of the triangulation 
Th- Further, let V/j and Wh be appropriate finite element subspaces of V and W satisfying 
LBB condition. For more examples of these spaces including Raviart-Thomas-Nedelec finite 
element spaces, Brezzi- Douglas-Marini spaces and Brezzi-Douglas-Fortin-Marini spaces, see 
HI]. 

The corresponding semidiscrete mixed finite element formulation is to seek a pair [uh-, (Th) '- 
(0,r] —)■ Wh X Vh such that 

(acT/i, Vh) - (u/,, V • V,,) = 0 V v/j G V/j, 

{uh,tu Wh) + (V • ah, Wh) = if, Wh) y Wh^Wh 


(3.1) 

(3.2) 
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with Uh{0) G Wh and Uh,tiO) G Wh to be defined later. Since Wh and V/j are finite 
dimensional, from (13.11) we can eliminate cTh in the discrete level by writing it in terms of 
Ufi- Therefore, substituting in (|3.2p . we obtain a second order linear system of ODEs and 
existence follows using ODE linear theory. Then, as a consequence of LBB condition and 
energy estimates, uniqueness can be proved easily and hence, we skip the proof. 

Set Bu = Uh — u and Scr = cTh — Erom (|2.4p - (l2.5l) and (I3.1l) - (|3.2p . and Co- satisfy 
the following equations 

(3.3) (aco-, v) - (e„, V • v) = ri(v) V v G V, 

(3.4) {eu,tt,w) + {V ■ ea,w) = r2{w) V tc G W, 

where the residuals ri and r 2 are given by 

ri(v) := (acT/,, v) - {uh,V ■ v), 


and 

r2{w) := {uh,tt,'w) + (V • ah,w) - if,w). 

Eollowing [2^, now introduce mixed elliptic reconstructions u{t) G Hq{Q) and cr(t) G V 
of Uh{t) and crh(t) for t G (0,T], respectively, as follows: for given Uh and let the mixed 
elliptic reconstructions u and a satisfy 

(3.5) (V • (ct - ah),w) = -r 2 {w), Vtc G W, 

(3.6) (a;(cr - cTft,), v) - (n - tt/i, V • v) =-ri(v), Vv G V. 

Using Theorem 4.3 (pp. 132) of [TO], one can verify that for a given Uh, crfi,ri and r 2 , the 
system (13.5h - (l3.6p has a unique pair of solution {u{t),a{t)} G TE x V, for t G (0,T]. Here 
elliptic reconstructions are assumed to be smooth in time. 

Note that ri(v;i) = 0 Vv/j G V/j, and r 2 (rc/i) = 0 Vwh G W^- Then, cr^ and are 
indeed mixed elliptic projections of a and u, respectively. 

Using mixed elliptic reconstructions, we now rewrite 

Bu := {u-u)-{u- Uh) =: iu - Vu, 

and 

:= (d- - cr) - (d- - (Th) =: - Vcr- 

An application of (I3.5l) - (|3.6p in (|3.3I) - (|3.4I1 yields 

(3.7) (a^^,v)-(C„,V-v) = OVvGV, 

(3.8) ■ C(T,U!) = {r]u,tt,w) \/w €W. 

With mixed elliptic reconstructions u and a satisfying (I3.5p - (l3.6h . apply (13.7p to check that 

(3.9) aa = —Vu. 
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Lemma 3.1. Let and ^a- satisfy (I3.7p - (l3.8p . Then, the following estimates hold: 

(3.10) Uu,tm + \\»^^"Ut)\\ < ||e.,t(0)|| + ||a'/'C.(0)||+2 r \\r^,,,uis)\\ds. 

Jo 

and 

||?«(t)|| < ||^n(0)|| + 2 / \\Vu,t{s)\\ ds. 

Jo 

Proof. Differentiate p3.7p with respect to t and set v = in the resulting equation to find 
that 

( 3 . 11 ) «,T,i,e^)-(eM>v-e^) = o. 

Choose w = (^u,t in P3.8p . Then, add the resulting equations to (13.lip to arrive at 

( 3 . 12 ) 1|(||^„_^||2 + ||« 1 / 2 ^^|| 2 ) ^ 

On integrating (I3.12p from 0 to t, a use of the Cauchy-Schwarz inequality yields 

||C„,i(t)f + < ||^n.i(0)f + ||a'/'C.(0)f + 2 T h„,«(s)|| ||C„,i(s)||ds. 

Jo 

Setting 

lll(«..<.{..)(<)lll = (l|a'''"{..(<)f + l|{.,i(*)P)'''". 

let t* G [0, t] be such that 

iii(e.,t,e.)(niii = max||i(e.,t,^.)(s)iii. 

0<s<t 

Then at time t = t*, equation ([3]) becomes 

|||(^n,,e.)(r)|||<|||(^M,e<.)(0)||| + 2 r \\VuAs)\\ds, 

Jo 

and hence, 

(3.13) |||(eM,e.)(t)lll < lll(e«.i,e.)(0)|||+2 r \\riu,u{s)\\ds. 

Jo 

This completes the proof of (j3.inp . Note that from (I3.13h . we obtain L°°(L^)-estimate of 
the displacement using ^(t) = ^(0) + /o^n,t(s) ds. Now in order to reduce the regularity, 
an integration tool which is a variant of Baker’s time testing procedure is used in a crucial 
way. To motivate our tool, integrate (13. 8p with respect to time to arrive at 

(3.14) w) + (V ■ ia-, w) = (^n,t(0), u;) + {riu,u w) - (?7„,t(0), u;), 

where fa- = fofa-is) ds. Choose rc = in (13.1411 and v = .^cr in (|3.7I) and adding the 
resulting equations to obtain 

14(||{„(t)|P + ||a‘''"j,(t)IP) = (e„,,(0), 
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Then, integrate with respect to time and use kick back arguments to arrive at 


||4(t)||<||U0)||+2 f\\vu,tmds. 

Jo 

This completes the rest of the proof. □ 

Assume that there exists a linear operator II/i : V ^ V/j such that V -Hh = Ph(V ■ ), 
where Ph '■ W ^ Wh is the L^-projection defined by 

(</) - Ph(t>, Wh) = ^ y Wh& Wh, (j>eW. 

Further, we assume that the finite element spaces satisfy the following properties: 


||v — II/ivll < C'//i^||v||r, 1 < r < + 1, llw' — Phw\\ < CiK^WwWr, 0 < r < ^ + 1. 

Note that for v € H(diu, Q) and w € L‘^{Q), the following properties hold true: 


{V ■ {v -Uh\'),Wh) = 0, Wh€Wh-, {w - PhW,V -Vh) = 0, Vh^Vh- 


Examples of spaces satisfying the above can be found in m- 

To prove the main theorem of this section, we need the following a posteriori estimates 
of r]u,Vu,t and rja- related to the mixed elliptic reconstructions (|3.5p - (l3.6l) . For a proof, see 

[El- 

Lemma 3.2. For Raviart-Thomas-Nedelec elements, there exists a positive constant C 
which depends only on the coefficient matrix A, the domain n, the shape regularity of the 
elements and polynomial degree i such that for £ = 0,1, 


(3.15) ||??«|| < C'(||/i^+^r 2 || + min \\h{acrh - VhWh)\\), 

and for j = 1,2, 


(3.16) 

and 

(3.17) 




dP 


< C 




19^2 


dP 


+ min 

wh&Wh 


h ( - '^hWh 


< C (^llhrall + \\ffi/‘^J{a(Th -1)110,+ \\h curlh {aa h)\\^ , 


where r 2 = {uh^tt “ / + V • cTh) is a residual and J{aah • t) denotes the jump of acrh ■ t 
across element edge E with t being the tangential unit vector along the edge E G F/j. 

Now, let £i{r 2 , a-h',Th), Ti(^^, 7^) and £ 2 {^ 2 ,(^h',Th) denote the terms on the 

right-hand sides of (I3.15p . (13.161) and (I3.17p . respectively. Then, using Lemmas 13.m3.21 we 
finally obtain the main theorem of this section as: 


























Theorem 3.1. Let {u,cr) be a solution of the mixed formulation (12.41) - (12. 5p and let {uh,crh) 
be a solution of the semidiscrete mixed formulation (I3II)-(I32]). Then the following a poste¬ 
riori estimates hold for £ = 0,1.' 


and 


I II I II 1/2 II 

< l|e«,t(0)|| ||a^/^e^(0)|| -\-£i{r2,t{0),crh,t{0);Th) 

+£2{r2{0),crh{0)]Th) + \\£i{r2,t, o-h,uTh)\\L^{o,T) 

+ 11^2(1-2, 7 ^)||x,oo(o,'r) + / £i{r2,tt,(Th,tt]Th) ds, 

Jo 


\^u\\L°°{o,T;L^{n)) ^ ||eM( 0 )||-I-£(i(r 2 ( 0 ), cr/j( 0 ); 7 ft) 


+ \\£i{r2,(Th;Th)\\L°°{o,T) + / £i{r2,t,crhTTh)ds. 

Jo 


4 Completely discrete scheme 

This section deals with a posteriori analysis for a completely discrete mixed approxima¬ 
tion based on backward differencing. 

Let 0 = tQ < ti < ... < t^ = T, In = {tn-i,tn\ and kn = tn — tn-i- For n G [0 : A^], let 
7)i be a refinement of a macro-triangulation which is a triangulation of the domain O that 
satishes the same conformity and shape regularity assumptions made on its refinements. 
Let 

hnix) := diam (K), where K ^ Tn and x € K, 

for all X € 17. Given two compatible triangulations Tn-i and Tn, i-e., they are rehnements 
of the same macro-triangulation, let Tn be the finest common coarsening of Tn and Tn-i, 
whose meshsize is given by hn '■= max(/in,see ([26], pp. 1655). 

We consider V)) and WJ^ defined over the triangulations T” as Raviart-Thomas finite 
element spaces of index ^ > 0 of Ti{div, £1) and respectively. Let : L‘^{£l) — WJ^ 

be the L^-projection dehned by 


{PJiw,clT) = {w,r) V 0 -GW,-. 

Given C/° = P^uo, hnd {([/”, E”)} with ([/", S’") G x V)) for n G [1 : iV] such that 

(4.1) w) + {V ■ S”, w) = ir, w) yw G WJi, 

(4.2) (aS",v)-([/",V-v) = 0 Vv G V)), 


where /*" is either chosen as point-wise value /”■(•) := f{tn,-) for n > 0 or the one through 
average value as in Remark 14.21 Here, the backward second and hrst hnite differences are 
given, respectively, by 




dtU^ 


- dtU^-^ 
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and 


dtU^ := 


jjn _ ^n-1 
kn 


Ph^l, 


for n = 1, • 
for n = 0. 


,N, 


Throughout the rest of the paper, we shall use the following notation: 

p^idtr) = 

Following m, we define for given a sequence of discrete values the time recon¬ 

struction y : [0, T] X n ^ iR or as 

(4.3) V{t) = V^ + {t- trr)dtV^ - < t < 

for n = 1, • • • ,N. Note that we have used the fact that dtV^ is well defined. 

We shall use the above C^-function V{t) such that for n = 0,1, • • • ,N, 

(4.4) V{t^) = V^, Vt{tn) = dtV^, Vu{t) = {l + ^l'^^)^^V^, 
for t € {tn-i,tn], where 

pt^(t) ;= 

Similarly, we define C^-functions U{t) and S(t) in time variable using the discrete 
sequences and {S"'}()Lq, respectively. 

As in Section 3, for given {U^, S"'}()Lq, we now define the mixed elliptic reconstructions 
u"’ € Hq (n) and d'"' G V at t = as: 


(4.5) 

(4.6) 


(V • - 51”), u;) = -r^(u;), w e W, 

(a(^- - E-), v) - (n- - V • v) = -r)^(v), v G V, 


where r”(v) := (q;S"',v) — ([/"■, V • v) and r^iw) := {PJ^{d^U"'),w) + (V • S^jtc) — {f'^,w). 

Since r”(v/i) = 0 Vv/^ G V)(, n > 0 and r'^{wh) = 0 Vtc/i G n > 1, in fact, 
'PP and [/” are mixed elliptic projections of d"” and iP at time t = tn, respectively. Now 
given {£("'}()Lq and {^^}n=ot we define the C^-functions u{t) and d'(t) in time t G (0, T], 
respectively, as 


(4.7) 
and 

(4.8) 


u{t) = iP + {t- tn)dtU^ - 


(t tn-l){tn t) „2~n 


di-lP, tn-l<t<tr, 


a{t) = + {t- tn)dcr'^ - —— — —8^0-^, tn-i <t<tn, 


provided that dtu^ and dfO-^ are well defined. 

For t G (0,T], the mixed elliptic reconstruction {n, d"} satisfies 

(V • (d- — 51), tc) = — r2(rc), w € W, 

{a{a - S), v) - (u - [/, V ■ v) = -ri(v), v G V, 


(4.9) 

(4.10) 
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where ri and r 2 are defined as C^-functions in time using ™ (14.3p . 

Again, set 

eu = {u-u) -{u- U) =: 

and 

= (ct - cr) - (or - E) =: ia- - r]„. 

Now, the pair {eu,ecr} satisfies 

(4.11) {eu,tt, w) + (V • e„,w) = {Utt,w) + (V • E, w) - (/, w). 

On splitting and Co-, we obtain from ()4.11l) 

(4.12) (e,,*,,u;) + (V-C^,u;) = ir^u,u,w) + {{I - P;:)Uu,w) + , Pj^w) 

+(V • (dr - dr*^), u;) + ir - f, w) Vw G W. 

Similarly, we also arrive at 

(4.13) (aC(T,v) - (4, V • v) = (Q;(di-- di-”),v) - (u - it”, V • v) Vv € V. 


Note that 

(4.14) CT - df” = (t - tn)dtd-^ + {k~^{tn - tf - (tn - t)^) O^a^, 
and 

(4.15) U-lP = [t- tn)dtlP + [k~^{tn - tf - [tn “ tf) . 

Now, it follows that 

(4.16) (a(df - dl-”),v) - (ft - u”, V • v) = (f - f„)|(a5tdl-"-,v) - (Stu", V • v)| 

+ {k~^{tn-t) - {tn - t)^) {(a5t^dr’",v) - V • v) |, 

and from (j4.5p with definition of ri(v), we find that 

(4.17) (adf",v)-(n”,V-v)=0 VvG V. 

From (I4.17p . the equation (|4.16l) takes the form 

(a(dr — df”), v) — (u — u”', V • v) = 0, 


and thus, (14.1311 becomes 


(4.18) 


(aC<T, v) - (^„, V • v) = 0 VvGV. 


Theorem 4.1. Let ( u ,( t ) and (C7, 51) be the solution of (|2.4l) - (l2.5p and (|4.II) - (14.21) . respec¬ 
tively. Then, the following estimates hold for t G {tn-i,tn] 


II^M(t)|| 


4 pt 

+ Wa^^^um < lien,i(0)|| + ||a'/2e.(0)|| +2Y,£iAt) + 2 / \\Vu,u{s)\\ ds, 

j=i 
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where 


£i,2it) 


^ 1 ptj pt 

J2 / \\iI-PDUu\\ds+ / \\{I-P^)Uu\\ds, 

j = 1 ^j — 1 in — 1 

n—1 

Y, \\^i^d^,W\\ds+ \\^P^jU^\\ds, 

T— 1 — 1 '^in — 1 


n-l /^2 


/c? 


f ||5t(r^2 - V • s^)|| + r^WdM - V ■ 5]- 


i=i 


^i,4(i) = 


+ / (i„-s)||9i(r^-V-5]’^)||ds 

'^in—1 

+ f ((t„ -»)" - ~ l|a?(rS - V ■ S")|| ds, 

in —1 V 'P' / 

ptj pt 

Y \\f-f\\ds+ \\r-f\\ds. 

j = \ — 1 in — 1 


Proof . Differentiate (14.18^ with respect to t. Then, choose r; = in the resulting equation 
and w = (,u,t in (14.121) to obtain for t E 

+ = iVu,U,^u,t) + {iI-P}:)Uu,fu,t) 


On integrating from 0 to t with t E we find that 

( 4 . 19 )^ ^ ds 

+ JUUt) + 

where 



".“1 ptj 

j=l 

f {{I-Pfi)Utuiu,t)ds, 

■^in—1 

dg2('^«,4) 

n —1 « 

■■= J2 1 h^id?U^,PlCu,t)ds+ / 

lP{d^,U^,Pfiiu,t)ds, 

ri —1 


■= Y [ - ^^)^^n,t)ds + 

j=l 

f {\/ ■ {a - a'^),^u,t)ds, 

^ ^n — l 


n—1 «+- 

_—, f^o _. j 

rt 


JUUt) -Y ^“ 4 ) + / ir - /, ^n,i) ds. 

j = \ ^j — l in — 1 
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Set 


El{t) := + 

and let at t = t* € (0, t] be such that 

Ei{t*) = max Ei{s). 

0<s<t 

Now, a use of the Cauchy-Schwarz inequality yields 

n-l 


cv'-yii ) 


E \\iI-PDUu\\ds+ [ \\{I - Pj:)Uu\\ds] E,{t*). 


and similarly, 

\Jl2i^u,t)\ < I E r \\f^^d^U^\\ds+ f \\^i'-^lU'-\\ds\ Er{P). 

\ j_l dtj-l Jtn-l I 

For J^g, we rewrite using (14.141) as 

(4.20) {V ■ {a - d-^),w) = {t - tn) ■ dta^, w) 

+ {k-\tn - - {tn - tf) (V • dja^^w) . 

From (3.6), we obtain 

(4.21) (V • ct", w) = -r^(n;) + (V • YP, w) , 
and therefore, for j = 1,2 

(4.22) {y-dia^,w) =-{div^){w) + {y-diYP^w) . 

On substituting (I4.22p for j = 1, 2 in Jj^g, we arrive at 


= E/ 

1=1 "'b-i 


,t)- 


- {kj\t^ - sf - (t,- - s)2) {(a2r^2)(Ci,,i) - (V ■ ) 


+ 


^n — 1 


(^{tn - s) {{dtr^){Ut) - (V • dtY^,Cu,t)} 

- {K\in - s)^ - {tn - sf) {{dfr^)i^u,t) - (V • 

Using the Cauchy-Schwarz inequality, it follows that 

/ f.2 ^3 \ 

|J[:3(Cm)| < E f dfri-V-Yf +:^ Ui(t*) 


i=i 


+ ( / {{tn-s)\\dt{T'^-V-Y^ 


^n — l 


+ {{tn - sf - k-\tn - s)^) \\d^{r^ - V • Y^)\\})Eft*). 
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For Jf 4 , we note that 


kr,4«..i)i = 


4-1 fti 

^ I 

j=i I'h- 


1’^ r 


- \ / I 

\i=i 


' {P-f,Cu,t)ds+ [ 

\\p-f\\ds+ f \\r-f\\ds]E,{n. 

On substituting the estimates of j = Ij •'' )4, in (I4.19p . we arrive at 

4 j-t 

Ei{t) < Ei{t*) < Ei{Q)+ 2'^Ei^j{t)+ 2 / \\r]u,tt{s)\\ ds. 

j=i 

This completes the rest of the proof. □ 

Remark 4.1. The term (r^ — V • S-^) can be replaeed by {d^U^ — f^). 

For obtaining L°°(L^) estimate for e^, we now integrate (j4.12p with respect to time from 
0 to t with t € {tn-i,tn], to arrive at 


(4.23) {^u,t,w) + \V ■ = {eu,ti0),w) + {r]u,t,w) 

-|-J£i(u)) -|- 

where 

n—1 


Jliiw) := {id-PDUu,w)ds+ f iiI-Pj:)Utt,w)ds, 

j — — — 1 

^ ^ /* ^3 C ^ 

J2,2{w) := X] / P{ciiU^,w)ds+ / fi'^{dfU'^,w)ds, 

j = \ —1 '^^n—1 

Pb —1 pfj pi 

J 2 pw) ■= '^ {V ■ {a-- a^),w)ds + {V ■ {a - a^),w)ds, 

j = \ ^j — 1 *^^n —1 

J£ 4 (u;) := ^ / {p-f,w)ds+ {r-f,w)ds. 

j = 1 '^^j — 1 '^^n—1 


Note that 


J2A = 


f {P-f,w)ds and J^pw) = f p{d'^U'^,w) ds 
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as ^ = 0. Further, since commutes with time derivative, can be written as 


n—1 




i=i 


+ ((/ - p;:)Ut{t),w) - ((/ - pj:)Ut{tn-i),w) 


n—l 


n—1 


(4.24) 


E ((' - -(i- pt')Ut(ti-i),w) + Y^iPl - pt'mti-i) 

i=i j=i 

+ ((/ - t/;) - ((/ - u;) 

n—l 

E - ((^ - ^°)O(0),«;) 

j=0 


Below, we prove one of the main results of this section. 

Theorem 4.2. Let {u,a) and (17,5]) be the solution of (|2.4l) - (l2.5p and (|4.II) - (14.21) , respec¬ 
tively. Then, fort G {tn-i,tn], the following estimate holds 

4 j-t 

||?«(t)|| < ||?„(0)||-F2^£'2j(t)-h2 / ||7/„,t(s)||cls, 

,=i 7o 

where £’ 2 , 1 , • • • ,£ 2,4 will be given in the following proof. 

Proof . Choose re = in (I4.23P and u = in (I4.18p . Then, add the resulting equations to 
obtain 


+ , fu) + {'nu,t,f,u) + J2,j{f,u)- 

i=i 

On integrating from 0 to t with t G (t„_i,t„], it follows that 

nt 4 

(4.25) ||^J|2 + ||«V2|^||2 = ||^,(0)||2 + 2(e„,t(0),e„) + 2 e.) ds + 2 ^ i7,-(t). 


i=i 


where 7£j(t) = / J 2 ,j{f,u)ds. Set 

Jo 

Elit) ■•= \\Ut)f + \\a^^HAt)\\\ 

and let t = t** € (0, t] be such that 


E 2 {t**) = max E 2 {s). 

0<s<t 
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Then, for tn-i < t <tn, we obtain from (14.2411 
n-i .tj i-i 

j=l -^b-i z=o 

rt n—1 

+1 - ^h)Ut{ti),Us)) ds - {{I-P^Mm,Us))ds 


'tn-l l^Q 

n-l ct 


+ E /' {iI-Pi)Ut{s),Us))ds+ [ {{I-Pj:)Utis),Us))ds. 

j = 1 — 1 tn—1 


Hence, 


\Kiit)\ < 


E E + (i - tn-iMPt^ - pim^\ 

j=l \ 1=0 

+ f' ||(/-P^)[/i(s)||ds] +(f-f„_i)||(pi-p°)9i[/°|| 


+fii(/ - + r ii(/ - p,")^t(5)iids 

in —1 

=: E2,im2{t**)- 

The second term can be written as 


E2{t* 


■n-l pt 

A,(«) . E/; 


j=l "'b-i \“'b-i 


+ 


H^{T){d^W,^u{s))dTj ds 
n"'{T){d^U^,Cu{s))dT] ds 


'in—1 V'iin—1 


for tn-l 


Since 


p 


/ fj.^{T)dT = -3kP 
'tj-i 

(s - tj-1/2) - -J 


we obtain 


l^2(t)| = 


^ ^ /*^j O ^ 1 ptj 

^E^7 / i^-''^j-i/2fi9tU^^^u{s))ds + -J2 kj{d^W,^u{s))d. 

j=i Pj-i j=i Pj-i 

t 


[ {s-tn-i/2f{d‘^U'^,^u{s))ds + jknf {8^ , Cu{s)) d. 

'^in —1 '^in —1 


< 


n—1 


Y.\\^PtU^\ + 

i=i 

=: £2,2{t)E2{t**). 


[{t - in- 1 / 2 )^ + ~ tn-l) + ) dfU 


kl 


E2{t**) 
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For the third term K^{t), we obtain for tn-i <t <tn and < s <tj, 

n—l j — l . 1 ^ 

Ksit) = '^ (V ■{a{T)-a^),^u{s))dTds 

j=i dtj-i dti_i 

pi ^ ^ 

+ (v ■ {a{T) - a^),Cu{s)) drds 

(V • (cr(T) - d-^),^u{s)) drds 
+ f f {V ■ {d-{T) - d-^),^u{s))dTds 

'^^n—1 '^^n—1 

= -fl + -^2 + -^3 + -^4- 

pi Ip 

Using (I4.20p and the fact that / — r)^ — (t; — r)^) dr = —we find that 

12 


n-i i-t 



j=l •'h-i 
ct rs 


"-1 rti 


k? 


n-1 rti i-1 


kf 


7=1 dtj-i 


^ Jf 12 

j=l -^h-l i=l 


lAl < 


n—l J*—1 / 7 2 / 3 

j=i 1=1 ^ 


E2it**). 


Similarly, 


\h\ < 


/jt2 . A:^ 

{t-tn-i)Y^ -^\\V ■ dta^W + -^\\V ■ d'^a^ 

j=i V 


For the I 3 and I 4 terms, we easily obtain 


and 


141 < 


V f(li^ _ S' 


j=i Ib-i 

n-l „i 


IIV • dta^W ds 


E2{t* 




< 


3 = 

■n-l /g 

3 


E2{t**) 


i=i 


S2(U 


|4| < {t-tn-l) 


yt|V-8ia-"|| + ^||V.8,V*|| E,(f 
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Now collect terms, replace V • cr-^ by —+ V • using (|4.2ip . and set £’ 2,3 := Mi + M 2 , 
where 


Ml = 

n-l 

E 

kj 

(E|ll«.(r 


i=i 

- 

\l = l 


+ {t 

tn 

-ofiiaK 

and 




Ms = 

n—l 

E 

1 

E|l|8?(r 


i=i 


\l = l 


K 

20 

(A - 

tn-l)\\dnr^2 

so that 





A;? /c? 


% 1 E f - V • soil I + 1 (A - An-l)f + f j WdM - V ■ S^' 


A:? kf 


1 E T^ll^'(r 2 - V • soil ) + ( (t - ^ ) WdKri - V ■ S^' 


|i^3(A)| <^2,3(A)i^2(A"). 

For the last term one can repeat previous arguments to arrive at 


l^4(t)| 


< 



£2A{t)E2{t**)- 


/(T)iidT + (t-t„_i) r iir 

in—1 


f{'r)\\dT 


E2{t**) 


On substituting in (I4.25p . it follows that 


||?u(A)|| < ||?u(0)|| + 2 ^£’2j(t) + 2 

i=i 



||f?«,t(s)|| ds, 


which completes the rest of the proof. □ 

In order to present the final theorem in this paper, we introduce some notations: For 
D := OT K, let 

£^{D) = ||ho(aSO + V;,17°)||i2p), 

£2iD) = (||C'r^||i2p) + ||h„(aS" + V,[/")||^2(^)), 

SSiD) = (||C5tr^||L2(^) + ||/i„5i(aS- + V,[/-)|U2p)), 

£2{D) = \\ho{adt'S^ + VhdtU^)\\mD), 

£l{D) = (||/iJ/V(aSO.t)||o,r„D + ||/iocurU(aSO)||i 2 p)), 

£^{D) = + l|AiyV(aS" • t)||o,r„D + IKcmh (aS-)||^ 2 (^)), 

S^D) = (||C5tr^||L2(^) + ||/i„5i(aS- + V,[/")||^2p)), 

SSiD) = (||C5MllL2(D) + l|/in52(«5]" + V,A/")||i2p)), 
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n - 

/ i'l 


and 




j=l 

n rU 


^1,2 = E/ Wd^dtU^Wds, 


j=i 




cn _ 
^1,4 — 


j=l \ 

j = l 


j -1 


- f{s)\\ds, 


fti 


cn _ 

^ 2,1 — 


E + / \\{i-pi)Ut{s)\\ds 

j=l \ 1=0 y 

+fc„ii(p,i - pji)dtu^+mi - pmtu% 


cn 

^ 2,2 


cn 

^2,3 


cn 

^2,4 


i=i 


kf. 


E % E Y Il'9t(r2 - V • S') + ^ 11^2(r' - V • S') 

j=i 1=1 ^ 

n .tj 

E^J / \\f-fi^)\\ds- 

7 = 1 ■^*7-1 


Using estimates of and r/o- in Theorems 14.11 and 14.21 the final theorem of this section can 

be written as 


Theorem 4.3. Let {u,a) be the solution of (|2.4p - (|2.5l) and (U, S) be the solution of (|4.1l) - 
Then for m € [1; -^]) the following estimates hold for the RT index £ = 0,1; 

m 4 

(4.26) \m-um\\ < \\eum\+Ci£^{n) + C2ST{n) + C3Y,kn£S + Y.P^U^^^ 

n=l i=l 

and 


||S™-<T(t^))|U-i < ||e„,i(0)|| + ||e^(0)|U-i+C4f°(fi) + C5f°(54) 

m 4 

(4.27) +Ce£^in) + C7£ip{n) + Cg 

n=l 2=1 

where Ci’s and Ci’s are constants which depend only on the eoeffieient matrix A, the domain 
Q, the shape regularity of the elements, polynomial degree i and interpolation constants. 

Remark 4.2. The last term in (|4.27|) . that is, £f^{£l) (also in (|4.26p . that is, £2 4 {£l)) 
measures the effect of approximating the forcing function f at discrete points in time. This 
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bound is of similar form as the bound in where the same diseretization has been 

used in the context of parabolic problems. However, a modification of /" as 


f 


1 



f{s)ds 


will specially improve the estimate in (|4.27p . 

Remark 4.3. The numerical implementation of the proposed a posteriori estimators in the 
adaptive algorithm deserves speeial attention and will be considered elsewhere. 


5 Conclusion 

The current work presents a first step towards true a posteriori estimate in the 
norm for mixed finite element approximations of second order wave equations. While Baker’s 
technique is usually used to derive estimates for the displacement u, in this pa¬ 

per, we resort to an application of integration for deriving these estimates. For adaptive 
algorithm, we need efficiency bounds, which would be an interesting direction for further 
research. Moreover, the numerical implementation of the adaptive algorithm based on the 
proposed estimators will be a part of our future work. 
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